omega = 1;
x = 0;
h = 0.08;
res1 = [trial1(omega,x,h),trial1(omega,x,h/2),trial1(omega,x,h/4)]
res2 = [trial1(omega,x,h),trial1(omega/2,x,h),trial1(omega/4,x,h)]
res3 = [trial2(omega,x,h),trial2(omega,x,h/2),trial2(omega,x,h/4)]
res4 = [trial2(omega,x,h),trial2(omega/2,x,h),trial2(omega/4,x,h)]
res5 = [trial3(omega,x,h)+1,trial3(omega,x,h/2)+1,trial3(omega,x,h/4)+1]
res6 = [trial3(omega,x,h)+omega,trial3(omega/2,x,h)+omega/8,trial3(omega/4,x,h)+omega/64]